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Motivation: Sign problem, Silver Blaze problem 
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Finite-density quantum chromodynamics (QCD) 


QCD 
Fundamental theory for quarks and gluons 


Neutron star 


e Cold and dense nuclear matter 


o 2Msun neutron star (2010) 


Neutron star merger 
o Gravitational-wave observations (2016~) (image from NASA) 


Reliable theoretical approach to equation of state must be developed! 


Zaco(T,p)= [DA Det(P(A, pa) +m) exp C Svu(). 


quark gluon 


Sign problem: Det(D(A, uqa) + m) Z 0 at uq 40. 
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Sign problem in finite-density QCD 


QCD & |QCD| 


ZQCD = [ra (det JD) e SYM | Z|QCD| = [ea [det P| eM, 


If these two were sufficiently similar, we have no practical problems. 
However, it was observed in lattice QCD simulation that at T = 0 


(e.g., Barbour et. al. (PRD 56 (1998) 7063)) 


nB QCD 


E > 
~My /2 -my /3 Pa 
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Baryon Silver Blaze problem 


The curious incident of the dog in the night-time (Holmes, Silver Blaze). 


Problem: Explain why ng = 0 for jug < my/3 using path integral. 
(Cohen, PRL 91 (2003) 222001) 


Current situation: For y < m; /2, the problem is solved. 
For lg > m4 /2, however, no one knows how to understand this. 
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Method: Path integral on Lefschetz thimbles 
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Sign problem of path integrals 


Consider the path integral: 


ja J «eec sta 


e S[x] is real = No sign problem. Monte Carlo works. 
e S|[x] is complex = Sign problem appears! 


If S[x] € C, eom S'[x] = 0 may have no real solutions z(t) € R. 


Idea: Complexify x(t) € C! 
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Lefschetz thimble for Airy integral 
Airy integral is given as 


3 
Ai(a) — EL (5 + 2 


Complexify the integration variable: z = x + iy. 


4 


Integrand on R, and on A 
(a — 1) 
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Rewrite the Airy integral 
There exists two Lefschetz thimbles J, (c = 1,2) for the Airy 


integral: 
Ai(a) Son f sae am 
= 5 — —+az). 
= J, 2T p 3 


ng: intersection number of the steepest ascent contour Ke and R. 
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Lefschetz decomposition formula 


Oscillatory integrals with many variables can be evaluated using the 
“steepest descent" cycles J,: (classical eom S'(2,) = 0) 


| dx e 560 = Y (Ko, y) d”z e80, 


o 


Jz are called Lefschetz thimbles, and Im[S] is constant on it: 


me =}, det) _ eS 


t——0o dt Oz? 


I= 20) 


(Ko, R}: intersection numbers of duals K, and R” 


(Ko = {2(0)|z(00) = zo y). 
[Pham, 1967; Kaminski, 1994; Witten, arXiv:1001.2933, 1009.6032] 
[Christoforetti et al. (PRD(2012)), Fujii et al. (JHEP 1310), etc.] 
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Previous studies: MC with one-thimble ansatz 
In previous studies, one has been used one-thimble ansatz: 


ll dex et) d^z e- 5), 
" Jo 


[Christoforetti et al. (PRD(2012)), Fujii et al. (JHEP 1310), etc.] 
Relativistic Bose gas: 


Thimble 2-(b) 5-5 " T 
3 C Langevin :-€-- > 
25 
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u (chemical potential in lattice unit) Fujii et al (JHEP 1310) 
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Analysis: Semi-classical analysis of the one-site Hubbard model 
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One-site Fermi Hubbard model 
One-site Hubbard model: 


A = Uff, — u(fu + f). 


Fock state gives the number density immediately: 


"d 


(YT, Hidaka, Hayata, 1509.07146)(cf. Monte Carlo with 1-thimble approx. gives a wrong result: 


Fujii, Kamata, Kikukawa,1509.08176, 1509.09141; Alexandru, Basar, Bedaque,1510.03258.) 
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Path integral for one-site model 
Effective Lagrangian of the one-site Hubbard model: 
p? 
L= oU + y* [0, — (U/2 + ip + u)] v. 
The path-integral expression is 


B ik d 2 052 
a=a — | ay (14 ef Gem)? 78/20 
27U R eee 


Fermion Det 


Integrand has complex phases causing the sign problem. 
ọ is an auxiliary field for the number density: 
(ñ) = Im(y)/U. 
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Sign problem and fermion determinant 


One-site Hubbard model: 
Uae —B(—U/2-1) oib)? 
Det |0, — jm enit =(l+e mete)”, 


Quark determinant in QCD: 


Det [aP , +m) — u] = N (A) lI (14-678 6;7u-i6)) (1 4 eles eH). 


ej>0 
where the spectrum of ya (1D, + m) is 
At) — €j (A) = id; (A) + (2n + Du. 


Minimal value of e(4) = m; /2. 
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Silver Blaze problem for y < —U/2, y < m4/2 


One-site Hubbard model: As BU >> 1 and —U/2 — u > 0, 


U . 
Det 2, = (» + 3 + ie) E (1 + e &CU/2- n ife y? e» 1. 


The sign problem almost disappears, so that J, ~ IR. 
Finite-density QCD: As 6 > oo and u < m4 /2, 


Det [ya(P4 +m) u] | II (1 + e-805-2)(1 + e- 802) 


= > [. 
Det[u(,-m) nao (+e) +eP) 


(Cohen PRL 91 (2003), Adams, PRD 70 (2004), Nagata et. al. PTEP 2012 ) 
The sign problem disappears by the reweighting method. 
=> Lefschetz thimbles ~ Original integration regions 
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Analysis 


Flows at /U < —0.5 (and u/U > 1/5) 


0.75 


0.25 


0 
Re(z) 


Figure: Flow at y/U = —1. Ja 


-0.25 


-0.5 


cR. 


; dz e 909, 


Number density: n, = 0 for u/U < —0.5, n, = 2 for u/U > 1.5. 


(YT, Hidaka, Hayata, 1509.07146) 


Z- 
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Silver Blaze problem for y > —U/2, y > m,/2 


One-site Hubbard model: At each real config., the magnitude is 
exponentially large: 


Det lo — (» + > + e) =0 (e8 U+#/2)) 


This large contributions must be canceled exactly in order for n = 0. 
Finite density QCD: The situation is almost the same, since 


Det (P(A, uj) + m) _ 
Det(D(A,0) +m) — 


I] ep8 qa), 


Re(AA) «uq 


but ng = 0 for y S my/3. 
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Analysis 


Flows at —0.5 < u/U < 1.5 


D 
[p 
Ss 
AAP E A SAR RRS 
ej 
sar 2S MAG GAMR 
2247: gd VM 
es, a 
207 £ G4 A RN NNR 
7 Ss, 
/ 
a a ES Nn 8 DR Je ai 
A IM ci 
E = 
AS ip ie RG ie ae onte 
^ 
y 
EEES => ==> 
ig n O a 
Be ae ANS So Sa S 
A M6 AAA A See eS 
eS o 
uu M. EA k JS 
id S 
com n a e a Oe 
AA a ee ee ee eT 
SS S 
~ 
E ee odis o abe ee ios Sr O 
== iui Eee uM Kcu MES ae” 
e e nou d mih e d tn k cc 
ie sia a ee) 
ss, A 
se o. A 7747 4— Te 
i 4 Lo 
—— s perenne oe eke ll 
= E 
UO SW de deo ae decir de gm m = 
SS p 
o OOS m cm cum 
` E eee 2 
"wo ox os we ccm 
eS À 
a Wa Ale umm mm n ele 
NC ao o 
Sos SNA cee wget le^ d 
NS ` bs 3 E 
NNN SO f m mE 
Mé---—^ 
NNN NAPE eC eer eee le 
e 
NON ON ON y TER EEE lA 
oco O O E e 
P. 
re) m a o 9 
M o pd S d 
S S 


Figure: Flow at u/U = 0 


wU p 1. 


This value is far away from n = Im (z)/U=0, 1, or 2. 


ra 


Complex saddle points lie on Im(z,,)/U 
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Curious incident of n in one-site Hubbard model 


We have a big difference bet. the exact result and naive expectation: 


05 0 0.5 1 15 ajo 


This is similar to what happens for QCD and |QCD|. 
p/U = -0.5 & py = Mx /2. 
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Complex classical solutions 


If GU > 1, the classical sol. 
for —0.5 < u/U < 1.5 
are labeled by m € Z: 


u[U x 1/2 


U 
Ea cR By + 27 mT. 


At these solutions, the classical actions become 


SE (B +3)’, 
Re (Sin — Sp) ~ 2m m?, 
Im S. c 2mm (5 4- 1). 
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Semiclassical partition function 


Using complex classical solutions Zm, let us calculate 


= 1 
Za = p» em 2 eme, (s (5 me 5) uud . 


This expression is valid for —1/2 S u/U $ 3/2. 
n4 + ny 


[U 
-1/2 0 1 a 
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Important interference among multiple thimbles 


Let us consider a “phase-quenched” multi-thimble approximation: 


Za = $e | = e 9 45(0, e 0727). 


o Lee-Yang zeros cannot appear at y/U — 0, 1. 
e One-thimble, or “phase-quenched”, result: n ~ u/U + 1/2. 


Consequence 


In order to describe the step functions, we need interference of 
complex phases among different Lefschetz thimbles. 


(cf. Particle Productions: Dumulu, Dunne, PRL 104 250402) 
(cf. Hidden Topological Angles: Behtash, Sulejmanpasic, Schafer, Ünsal, PRL 115 041601) 
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Numerical results 


Results for GU = 30: (1, 3, 5-thimble approx.: Jo, Jo U Ji, and Jo U J41 U Ja ) 
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Necessary number of Lefschetz thimbles ~ 8U/(27). 


(YT, Hidaka, Hayata, NJP 18 (2016) 033002, arXiv:1509.07146[hep-th]) 
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Review 


Review: Practical algorithm for simulating multiple Lefschetz thimbles 


Yuya Tanizaki (RBRC) 


Lefschetz-thimble path integral 


DAG 


Review 


Possible concerns for practical applications 


In this study, we emphasized the importance of interference among 
Lefschetz thimbles. 


Then, the following becomes important: 
e Find all contributing complex saddle points, 
e Construct Lefschetz thimbles around those saddle points, 
e Evaluate integration on each Lefschetz thimbles, and 
e Sum up those results. 
We need some machinery to do them automatically. 
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Review 


Neat idea for multiple thimble simulation 


"Y . dz os 
Deform the original cycle R” by the gradient flow, Ti (5) 
Z 
(Alexandru, Basar, Bedaque, Ridgway, Warrington, JHEP (2016)) 
Im(z) 
Re(z) 
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Review 


Neat idea for multiple thimble simulation 


d oS 
Deform the original cycle IR” by the gradient flow, E == (5) 


(Alexandru, Basar, Bedaque, Ridgway, Warrington, JHEP (2016)) 


Im(z) 
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Review 


Neat idea for multiple thimble simulation 


Oz 


i d os 
Deform the original cycle R” by the gradient flow, E == ( ) 


(Alexandru, Basar, Bedaque, Ridgway, Warrington, JHEP (2016)) 


large action 
X small action varying Im[S] 
A -constant Im[S] ~constantAm[S] 


Thanks to Gokce for nice figure! 
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Review 


Formulation 
Let us fix a flow time T', and define 


UP em far x) € C” aha) = (5) z(0;::z) =x €] i 


A 


By construction, z(T;-) : IR^ 2» J(T) and 


] seo = i dz e-50) 
B I(T) 


i d"z det (22) e- S62). 
n Ox) 


— One can do usual Monte Carlo + reweighting by regarding 


Sería) := S(z(T;x)) — In L (2) 


Ox) 


as the effective classical action. 
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Review 


Real-time dynamics 
This method is applied to Schwinger-Keldysh path integral for 


u l l3 4 
s- fa (s go 2 


03? a LEN ur Re[Gr (t.t)] 
= EN — Im[Grx(tt')] é $ 


0.0 0.5 1.0 1.5 2.0 


Feynman propagators at 8 = 0.8. Tgow = 0.2. (Alexandru,Basar, Bedaque, Vartak, 


Warrington, arXiv:1605.08040) 
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Summary & Conclusion 


Summary and conclusion 
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Summary & Conclusion 


Summary and Conclusion 


e Lefschetz-thimble method gives a suitable framework for 
saddle-point analysis even if S(@) takes complex values. 


e One-site Hubbard model is a nice toy model to play with the 
sign problem. 


e Interference of complex phases among Lefschetz thimbles plays a 
pivotal role for the (baryon) Silver Blaze problem. 


e Recent developments may enable us to study nonperturbative 
field theories with the sign problem. 
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